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Abstract 

We note that Fujikawa's proposal of generalization of the Ginsparg- 
Wilson relation is equivalent to setting R = {a^D) 2k in the original 
Ginsparg- Wilson relation D75 + 75D = 2aDR / y^D. An explicit realiza- 
tion of D follows from the Overlap construction. The general properties 
of D are derived. The chiral properties of these higher-order ( k > 
) realizations of Overlap Dirac operator are compared to those of the 
Neuberger-Dirac operator ( k = ), in terms of the fermion propagator, 
the axial anomaly and the fermion determinant in a background gauge 
field. Our present results ( up to lattice size 16 x 16 ) indicate that the 
chiral properties of the Neuberger-Dirac operator are better than those 
of higher-order ones. 



PACS numbers: 11. 15. Ha ( Lattice gauge theory ), 11.30.Rd ( Chiral 
symmetry ) 

Keywords : Lattice gauge theory, Chiral fermion, Ginsparg- Wilson re- 
lation, Axial anomaly. 



1 Introduction 



During the last two years, it has become clear that the proper way to formulate 
chiral fermions on the lattice is to impose the exact chiral symmetry on the 
lattice, namely, the Ginsparg- Wilson relation |1| 

D l5 + l5 D = 2aDR l5 D , (1) 

where D is the lattice Dirac operator, a is the lattice spacing, and R is a 
positive definite Hermitian operator which commutes with 75. Equation ([!]) 
should be regarded as a generalized chiral symmetry which contains the usual 
chiral symmetry in the continuum limit ( a — *■ ). However, it should be noted 
that this exact chiral symmetry ( R — 1 ) had been existing in the Overlap 
formalism |2|, ||, even before the GW relation was rediscovered. Therefore, 
unless one can explicitly construct a GW Dirac operator D without using 
the Overlap, and such a D satisfies all physical requirements; otherwise it is 
unlikely that the GW relation would turn out to be more fundamental than 
the Overlap. 

Recently, Fujikawa @] proposed a generalization of the Ginsparg- Wilson 
relation as 

75(75^) + (75^)75 = 2a 2k+ \ 1 ,Df k+ \ k = 0,1,2,- ■■ (2) 

Multiplying both sides of (0) by 75, we obtain 

75 £ + £75 = 2aD(a l5 D) 2 %D, k = 0, 1, 2, • • • (3) 

which is equivalent to the original GW relation ([!]) with 

R = (a l5 D) 2k , fc = 0,1,2,-.. (4) 

It can be shown that R is Hermitian and commutes with 75 for any D which is 
75-Hermitian and satisfies (|3|) [ the proof is below Eq. fliTf ) ] . The motivation 
of considering k > in (f|) is to improve the chiral symmetry at small lattice 
spacings. Further, Fujikawa has constructed a sequence ( k > ) of these GW 
Dirac operators based on the Neuberger-Dirac operator ( k — ) ||. However, 
the price one has to pay for the improved chiral symmetry is a less localized D, 
since in the limit k — > 00, D must tend to a chirally symmetric and nonlocal 
D c , as a consequence of the Nielson-Ninomiya no-go theorem ||. Therefore, it 
is not clear whether one may have any advantages in practice by considering 
k > in (|). Nevertheless, from a theoretical viewpoint, it is interesting to see 
how one can construct a sequence ( k — 1, 2, • - ■ ) of topologically proper D 
satisfying the GW relation (|3|), in addition to the Neuberger-Dirac operator ( 
k = 0). 

In this paper, we examine several aspects of Fujikawa's proposal. In section 
2, we derive the analytical properties of the GW Dirac operator satisfying ([[]) 



1 



with R = {a^j^D) 2k . Then, in section 3, we analyze the construction of higher- 
order Overlap Dirac operators, and derive their general properties. In section 
4, we compare the chiral properties of the higher-order ( k — 1, 2 ) Overlap 
Dirac operators to those of the Neuberger-Dirac operator, by computing the 
fermion propagator, the axial anomaly and the fermion determinant in two- 
dimensional background U(l) gauge fields. Finally, we discuss and conclude 
in section 5. 



2 General analytical properties 

In this section, we begin with general considerations of the GW relation, and 
then derive the analytical properties of the GW Dirac operators satisfying ([]]) 
with R = (a^D) 2k . 

In general, one can assume that the lattice Dirac operator D satisfies the 
Ginsparg- Wilson relation in the form 

D l5 f(D) + g(D) l5 D = , (5) 

where / and g are any analytic functions. Then the fermionic action Af = 
ipDip is invariant under the global chiral transformation 

^ _> exp[e l5 f(D)) i> (6) 
tP - $ exp[0<7(£>)7 5 ] (7) 

where 9 is a global parameter. 

If we set f(D) = 1 — aRD and g{D) = 1 — aDR, then (JD becomes 

D lb + l5 D = aD(R l5 + l5 R)D (8) 

where R is any operator. 

Since the massless Dirac operator in continuum is chirally symmetric ( 
T>7 5 + 75D = ) and antihermitian ( T)> = —T> ), so it is 75-Hermitian 
( J)' = 75l , 75 ). Thus, we require that the lattice Dirac operator D also 
preserves this symmetry at any lattice spacing, i.e., 

= 75^75 ■ (9) 

Then multiplying (^) by 75 and using (^), we obtain 

D ] + D = aD\R + 75^75)^ = aD(R + l5 R lb )D^ . (10) 

Evidently, only the part of R which commutes with 75 can enter (0). Recall 
that any operator R can be decomposed into two parts as 

R = \{R + 75^75) + \{R- 75^75) 
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where the first ( second ) term on r.h.s. commutes ( anticommutes ) with 75. 
Therefore, without loss, one can assume that R commutes with 75. Thus, ( |10D 
becomes 

D^ + D = 2aD^RD (11) 

Taking the adjoint of (|TTD , we immediately obtain 

R = R ] . (12) 

Then (§) becomes the usual GW relation 

D l5 + l5 D = 2aDR l5 D , (13) 

where R is any Hermitian operator which commutes with 75. 
Fujikawa's proposal Q is equivalent to setting 

R=(a l5 D) 2k , fc = 0,1,2,-.. (14) 

in the usual GW relation (|T3|). Then fliTf ) becomes 

^75 + 75^0 = 2aD(a l5 D) 2k l5 D . (15) 

It is obvious that R is Hermitian since D is 75-Hermitian. Note that 

l5 (a l5 D) 2 = (a 75J D) 2 75 (16) 

since 

75 (cr^D) 2 

= a 2 75 (7 5 £ + ^7s) (IbD) - a 2 (j 5 D)j 5 (j 5 D + D l5 ) + a 2 ( l5 D) l5 (D l5 ) 
= 2(a l5 D) 2k+2 - 2(a l5 D) 2k+2 + (a l5 D) 2 l5 
= (a7 5 D) 2 7 5 

where (|l^) has been used in the second equality. Then it follows that R = 
(a75-D) 2fc commutes with 75. 
Further, Eq. (|16D gives 

l 5 (l 5 D) 2 l5 = ( l5 D) 2 (17) 

which yields 

DD ] = D^D, D is normal (18) 

where @ has been used. Since D is normal, D and D> have common eigen- 
functions and their eigenvalues are either real or come in complex conjugate 
pairs, 

D<P S = A S S , (19) 
D^ 8 = A*0 S , (20) 
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(24) 



and the eigenf unctions {<f) s } form a complete orthonormal set. Then the gen- 
eral analytical properties of the eigenmodes as derived in Ref. || ( Eqs. (37), 
(38) and (41) in Ref. § ) for R = 1/2 also hold for the R in (|TJ). For the 
sake of completeness, we outline the derivations as follows. Writing R as 

R = (a l5 D) 2k = a 2k ( l5 D l5 D) k = a 2k (D^D) k , (21) 

and applying Eq. (|IID to <p s , we obtain the eigenvalue equation 

X s + A: = 2a 2fe+1 (A:) fc+1 (A s ) fc+1 . (22) 

Multiplying both sides of (0) by (X s ) k {X*) k , we get 

(X s ) k+1 (X* s ) k + (X s ) k (Xl) k+1 = 2a 2k+1 (X* s ) 2k+1 (X s ) 2k+1 (23) 

which can be rewritten as 

1 

2 

Thus the eigenvalues in the form z = (aX s ) k+1 (aX*) k fall on the circle centered 
at 1/2 with radius 1/2. The real eigenvalues ( if any ) of D are at A = ( z = 
) and A = a -1 ( z = 1 ). Writing aX s = rexp(i8), we have z = r 2k+1 exp(i6) . 
Thus, for k = 0, the eigenvalues of D fall on the circle centered at 1/2 with 
radius 1/2, while for k > 0, on the deformed circle stretched symmetrically 
in ±y directions with fixed points at zero and a" 1 , bounded inside the region 
< x < a -1 . 

From (|S|) and ([ZUD, we obtain 

Di 5 <J) s = A: 75 s • (25) 
Multiplying both sides of ( p5|) by (p\ and using the adjoint of (|20|), we get 

A^750 s = A>t 7 50 s (26) 

This implies that the chirality of any complex eigenmode is zero, 

Xs = <t>l l5 <t> s = if A s ^ A:. (27) 

If X s is real ( zero or a -1 ), then Eqs. ( p5|) and ( |1~9D imply that <p s has definite 
chirality +1 or — 1 : 

750 s = ±<f>„ if A s = A:. (28) 

A useful property of chirality is that the total chirality of all eigenmodes must 
vanish, 



E£££[tf(*)r7ftf(*) 

s x a p 

££75^ = (29) 



P 
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where the completeness relation 

££[C(z)]*0f(*) = ^ (so) 

X s 

has been used. Since the chirality of any complex eigenmode is zero, then Eq. 
( p9|) gives the chirality sum rule || for real eigenmodes, 

iV + + n + = JV_ + n_ (31) 



where n + (nS) denotes the number of zero modes of positive ( negative ) chi- 
rality, and N + (NS) the number of nonzero real ( a -1 ) eigenmodes of positive 
( negative ) chirality. Then we immediately see that any zero mode must be 
accompanied by a real ( a~ l ) eigenmode with opposite chirality, and the index 
of D is 



index(L)) = n_ - n + = -(AL - iV + ) . (32) 

It should be emphasized that the chiral properties (p8l), and the chirality 
sum rule ( |3~1"D hold for any normal D satisfying the 75-Hermiticity, as shown 
in Ref. ||. However, in nontrivial gauge backgrounds, whether D possesses 
any zero modes or not relies on the topological characteristics [J7| of D, which 
cannot be guaranteed by the conditions such as the locality, free of species 
doublings, correct continuum behavior, 75-Hermiticity and the GW relation. 



3 Higher-order realization of Overlap Dirac 
operator 

Now the task is to construct a topologically proper D which is local, free 
of species doubling, 75-Hermitian, having correct continuum behavior, and 
satisfies the GW relation (|3[). So far, the only viable way to construct a 
topologically proper D is the Overlap, 

£>= -1(1 + 756), e 2 = I (33) 
2a 

which satisfies the GW relation ([!]) with R = 1. There are many different 
ways to implement the Hermitian e in fl3"3]). However, it is required to be able 



to capture the topology of the gauge background. That means, one-half of 
the difference of the numbers ( h± ) of positive ( +1 ) and negative ( — 1 ) 
eigenvalues of e is equal to the background topological charge Q, 

J2tr(a l5 D(x,x)) = \ y £ti{e) = \{h + - h.) = Q , (34) 
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where tr denotes the trace over the Dirac and color space. Otherwise, the axial 
anomaly of D cannot agree with the topological charge density in a nontrivial 
gauge background. Henceforth, we shall regard any lattice Dirac operator 
which is constructed through the general form of Overlap (|33|) as a realization 
of the Overlap Dirac operator. 

An explicit realization of e in fl33|) is the Neuberger-Dirac operator || with 

H,r (35) 



where 

H w = 75 (Av - ^ocT 1 ), < m < 2r w , (36) 
Dw = liitu + W, Dw '■ massless Wilson-Dirac operator , (37) 

t»(x,y) = ^ [U„(x)S x+ a, y - Ul{y)S x ^ y ] , (38) 
W (?>V) = r ^T,[ 2 ^v- U ^)^v- U l(y)^-^ ■ ( 39 ) 
( The Wilson parameter r w is usually set to one. ) 

7 " = ( °J o ) • 

and 

+ o v o\ = 25^ . 

Note that the parameter m plays a crucial role in detecting the topology of 
the gauge background. 

Now the problem is how to generalize this construction to D satisfying the 
GW relation with R = (a r )z ) D) 2k for k > 0. We can multiply both sides of the 
GW relation ([!]) by R and redefine D' = RD, then we have 

£'7s + 75^' = 2D' l5 D' (40) 

where D' = RD is 75-hermitian since 

D't = D*R = 75^75^ = 75^^75 = #75^75 = J 5 RDj 5 = j 5 D' l5 . 



Now ( f40D is in the same form of the GW relation with R — 1. Thus, one can 
construct D' in the same way as the Overlap 

D' = RD = ^-(l + l5 e), e 2 = l, (41) 
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provided that a proper realization of e can be obtained. An explicit construc- 
tion based on the Neuberger-Dirac operator ( k — ) has been generalized to 
higher-orders ( k > ) by Fujikawa [|]]. 

In the following, we formulate Fujikawa's construction in a more transpar- 
ent way. Using R^ 5 = 75-R, we can write 

D' = RD = (a l5 D) 2k D = (a l5 D) 2k l5l5 D = a - l l5 (a lb D) 2k+l , (42) 

which yields 

D = a- 1 75 (a 75J D') 1/{2fc+1) , (43) 

where the (2k + l)-th real root of the Hermitian operator a^^D' is assumed. 

Then (|43| ) suggests that if the e in ([|l]) is expressed in terms of a Hermitian 
operator H, 

H (44) 



VH 2 

then H is required to be proportional to (7 ^T> fJi ) 2k+1 plus higher-order terms in 
the continuum limit such that D behaves as after taking the (2k + l)-th 
root in (f43|). Thus, H must contain the term (•y fl t fl ) 2k+1 , where 7 M t M is the 
naive lattice fermion operator defined in (p8|). Then additional terms must 
be required in order to remove the species doublers in the term (7 A1 t At ) 2A:+1 . 
So, we add the Wilson term to the (2k + l)-th power, i.e., W 2k+l . Finally, a 
negative mass term — (moa _1 ) 2fc+1 is inserted such that e is able to detect the 
topological charge Q of the gauge background, i.e., 

kE*(-7m)=Q- ( 45 ) 



2^ Wh\ 

Putting all these terms together, we have 

H = 75 [(7^) 2fc+1 + W 2k+1 - (m,a- l ) 2k+1 ] , (46) 
which, at k — 0, reduces to the H w in Eq. ([361). Then ([4~3|) can be rewritten as 



• lx l/(2fc+l) / tt \1/(2A:+1) 

d (-) 75 h + 7ff) ' ^ 



where H is defined in fl4"6p. This is the higher-order realization of Overlap 
Dirac operator, as constructed by Fujikawa 0]. 

Next we derive some general properties of the Overlap Dirac operator D 
defined in (g7|). 

The fermion propagator Sf(x,u) is defined by 

S F (x,y) = i| Hd${z)dif>{z) e"^ ^(a;)^(y) 
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where 



Z= Y[dip(z)dip(z)e 



■ipDip 



In a background gauge field of zero topological charge ( Q = ), the fermion 
propagator is 



S F (x,y) = D 1 (x,y) . 



(4£ 



In the naive continuum limit with r w a~ l and mod -1 kept finite, the ( free ) 
fermion propagator in momentum space can be obtained after some straight- 
forward algebras, 



S F (p) = 075 



75 1 + 



1 



(7A) 2fc+1 



Tip) 



l/(2fe+l) 



where 



ia 1 sin(p M a) 
a' 2 sin 2 (p i 



W(p) = r w a 1 [1 - cos(p M a)] 



u(p) 
iV(p) 

r(p) 



y/(t 2 ) 2k+1 + u 2 (p) , 
iV(p) — iz(p) . 



(49) 



(50) 
(51) 

(52) 

(53) 
(54) 
(55) 



Around p ~ 0, £ M ~ T(p) ~ 2{rriQa 1 ) 2fc+1 ) and the fermion propagator 
becomes 



&(p) 



«75 



75 1 + 



1 



-2(moa 



l/(2fe+l) 



(W/J 2fc+1 
2 i/( 2fc+ i) mo _JL_ + a + a (l _ 5 M )E fe (ap) 



(56) 
(57) 



where except for k = 0, a momentum-dependent term denoted by aEfc(ap) is 
present in the scalar part, which may lead to additive mass renormalization. 
Evidently, we have to fix the value of mo to 



m 



j^i/pfc+i) 
2, 



(58) 



such that in the limit ( a — > ) the fermion propagator fl5"7|) agrees with the 
continuum propagator. 
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For itlq e (0, 2r w ), on a <i-dimensional lattice ( d = even ), at any one the 
2 d — 1 corners of the Brillouin zone [ i.e., ap = (it, 0, • • • , 0), (0, ir, ■ ■ ■ , 0), • • •, 
(71", 7r, • • ■ , 7r) ], we have AT(p) = w(p) > 0, thus T(p) = — u{p) = 0, and all 
doubled modes are decoupled from the fermion propagator (|49| ) . 

In general, we consider an arbitrary value of itlq. At the origin ( p = ) 
and the 2 d — 1 corners of the Brillouin zone, we have sin(p M a) = 0, so T(p) 
becomes 

T(p) = \u(p)\-u(p) , (59) 

where the possible values of u(p) are : 

u{p) = -(moa" 1 ) 2 ^ 1 , 

(2r^a- 1 ) 2fc+1 - (moa- 1 ) 2 ^ 1 , 



(2rfr w a~ 1 ) 2fc+1 - (moa' 1 ) 2 ^ 1 . 

Here the first value of u{p) corresponds to all components of p equal to zero, 
the second value to one of components equal to ir/a, and so on, and the last 
value to all components equal to n/a. Note that 

{2nr w a- l ) 2k+1 - {m a~ l ) 2k+l 
= {2nr w a~ l — m a^ 1 )[(2nr w a^ 1 ) 2k H h (m a _1 ) 2fc ], n = 0, 1, • • • , d. 

Therefore the sign of u(p) is independent of the order k. From Eq. (|59D, we see 
that if u(p) > 0, then T{p) = 0, and this doubled mode is decoupled from the 
fermion propagator fl49|) for any order k. Since the chiral charge of a doubled 
mode is equal to (— l) n , where n is the number of momentum components 
equal to n/a, then the total chiral charge Qa of all massless ( primary and 
doubled ) fermion modes contributing to the fermion propagator ( f£9"D can be 
determined. Then index(D) = Q$Q for a gauge background with topological 
charge Q. Thus, in the naive continuum limit, the index of D ( as a function of 
mo ) is independent of the order k, same as the index of the Neuberger-Dirac 
operator ( k = ), which has been determined in Ref. ||. Explicitly, 

index[D(m<,)] = ^ for n = 1, • • • , d ; (60) 

0, otherwise. 
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In particular, for d = 4, 



index [D(m )] = < 



Q, 


< 


m < 2r w , 


-3Q, 




< m < 4r w 


3Q, 




< m < 6r w 


-Q, 


6r„, 


< m < 8r w 


0, 


otherwise. 



(61) 



For the Neuberger-Dirac operator ( k — ), there exists an exact discrete 
symmetry of the index on any finite lattice with even number of sites in each 
dimension || 

index[D(m )] = — index [D{2dr w — m )], for k = 0, (62) 

which holds for any background gauge configuration. However, for higher- 
order Overlap Dirac operators ( k > ), this discrete symmetry is not exact 
on a finite lattice; only in the naive continuum limit, this discrete symmetry 
can be realized as in (|60|). Nevertheless, at m can be shown that 

the index is exactly zero for all k, 

index[_D(m = dr w )] = 0, for all k > , (63) 

on any finite lattice with even number of sites in each dimension, and for any 
background gauge configuration. 



4 Tests 

In this section, we compare the chiral properties of the higher-order Overlap 
Dirac operators ( k = 1, 2 ) to those of the Neuberger-Dirac operator ( m = 1 
and R — 1/2 ), by computing the fermion propagator, the axial anomaly and 
the fermion determinant in two-dimensional background U(l) gauge fields. 
Our notations for the two-dimensional background gauge field are the same as 
those in Ref. © ( Eqs. (7)-(ll) in Ref. |f ). 

Note that the Neuberger-Dirac operator is conventionally written as 

D = a- 1 (l + V) = a- l (l + 1 r -^L), 77io = l; (64) 




which satisfies the GW relation with R = 1/2. However, the zeroth order 
k = ) Overlap Dirac operator is 
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which satisfies the GW relation with R = 1. In the following, we shall use the 
Neuberger-Dirac operator (|64D in place of the zeroth order Overlap d65|) . All 
numerical results for the k = case are obtained using the Neuberger-Dirac 
operator (|B3D rather than fl65p . 

First of all, we checked that the eigenvalues of a higher-order ( k = 1, 2 ) 
Overlap Dirac operator fall on the deformed circle which is stretched symmet- 
rically in ±y directions. In a nontrivial gauge background, the real eigenmodes 
( zero and a -1 ) have definite chirality and satisfy the chirality sum rule fl31]), 
and each complex eigenmode has zero chirality fl27|) . The Atiyah-Singer index 
theorem ( n_ — n + = Q ) is satisfied in all cases ( k — 0, 1, 2 ) for gauge 
configurations fulfiling the topological bound |7] 

a 2 |p(x)| < e l ~ 0.28 Vx (66) 

where a 2 p(x) is the topological charge inside the unit square of area a 2 centered 
at x, 

1 rXx+a/2 rX2+a/2 \ 

p( x ) = ~2 d vi d vi tt f ^v) ■ ( 67 ) 

a z Jxi-a/2 Jx 2 -a/2 2,71 

The value of mo in the higher-order ( k > ) Overlap Dirac operator is fixed 
according to Eq. (0), niQ = 2~ 1 /( 2k+1 \ while mo = 1 for the Neuberger-Dirac 
operator. 



4.1 Fermion propagator 

In the following, we first compute the free fermion propagator on a two dimen- 
sional lattice, for the Neuberger-Dirac operator and higher-order ( k = 1, 2 
) Overlap Dirac operators respectively. We compare them to the exact so- 
lution of the massless fermion propagator on the torus. Then we turn on a 
background gauge field to examine the behaviors of the scalar part So and the 
pseudoscalar part S 5 in the higher order ( k — 1, 2 ) fermion propagators. 
In general, the free fermion propagator can be written as 

S F (x) = S (x) + j^S^x) , (68) 

where So(x) = for the massless fermion in continuum; and Sq(x) = (a/2)5 Xt o 
for the Neuberger-Dirac operator. 



First, we examine the S^(x) components in (|68|). In Table 1, we list the 
component Si(x) along the diagonal ( x\ = £2 ) of a 16 x 16 lattice with 
antiperiodic boundary conditions. One of the end points of the propagator 
is fixed at the origin, while the other end point is located at a site along the 
diagonal. Note that, by symmetry, 5*2(0;) = Si(x) along the diagonal. From 
the data in Table 1, we immediately see that in all cases ( Neuberger, k = 1, 2 
), Si(x) agrees very well with the exact solution on the torus. In general, all 
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x i 


Exact 


Neuberger 


k = 1 


k = 2 


0.0 


0.0 


OOOO 


0000 


0000 


0000 


1.0 


1.0 


0.0796 


0689 


0.0894 


0885 


9 


2 


0396 


0374 


0423 


0449 


3.0 


3.0 


0259 


0.0249 


0226 


0.0212 


4 


4 


n 01 84 


01 79 


01 77 


01 79 


5.0 


5.0 


0.0131 


0.0127 


0.0128 


0.0126 


6 


fi 


8 59 x 1 0~ 3 


8 32 x 1 -3 


8 32 x 1 0~ 3 


8 38 x 1 0~ 3 

O .UO -A. 1U 


7.0 


7.0 


4 28 x 10~ 3 


4.13 x 10~ 3 


4.10 x 10 -3 


4.17 x 10 -3 


8.0 


8.0 


0.0000 


0.0000 


0.0000 


0.0000 


9.0 


9.0 


-4.28 x 10~ 3 


-4.13 x 10~ 3 


-4.10 x 10~ 3 


-4.17 x 10~ 3 


10.0 


10.0 


-8.59 x 10" 3 


-8.32 x 10" 3 


-8.32 x 10" 3 


-8.38 x 10" 3 


11.0 


11.0 


-0.0131 


-0.0127 


-0.0128 


-0.0126 


12.0 


12.0 


-0.0184 


-0.0179 


-0.0177 


-0.0179 


13.0 


13.0 


-0.0259 


-0.0249 


-0.0226 


-0.0212 


14.0 


14.0 


-0.0396 


-0.0374 


-0.0423 


-0.0449 


15.0 


15.0 


-0.0796 


-0.0689 


-0.0894 


-0.0885 



Table 1: The free fermion propagators on a 16 x 16 lattice with antiperiodic 
boundary conditions. One of the end points of the propagator is fixed at the 
origin, while the other end point is at one of the sites along the diagonal ( 
xi = x 2 ). 

S/j, components of the free fermion propagators are in good agreement with the 
exact solution for any x = (xi,^). 

Next, we examine the scalar part Sq(x) in the free fermion propagator of 
the higher-order ( k = 1,2 ) Overlap Dirac operators (|47|). In Table 2, we 
list Sq(x) along the diagonal of the 16 x 16 lattice with antiperiodic boundary 
conditions. From the data in Table 2, we see that Sq(x) is local for both k — 1 
and k = 2. However, we note that Sq(x) in the second order ( k = 2 ) case 
is less localized than that of the first order ( k — 1 ), as expected. It seems 
that |So(x)| in both ( k — 1, 2 ) orders can be fitted by an exponentially decay 
function for < \x\ < 8y2. 

In a background gauge field, the fermion propagator can be written as 

S (x,y) + S 5 (x,y) S R (x,y) \ , gg , 

S L (x,y) S (x,y) - S 5 (x,y) J ' 

where So(x,y) = S^(x,y) = for the massless fermion in continuum; and 
So(x,y) = (a/2)5 Xt y and S$(x,y) = for the Neuberger- Dirac operator. How- 
ever, for higher-order ( k > ) Overlap Dirac operators, both So(x,y) and 
S 5 (x,y) are not proportional to S x>y . If S (x,y) ( S 5 (x,y) ) turns out to be 
nonlocal, then it would cause additive mass renormalization and the poles in 
the fermion propagator will be shifted accordingly. 



S F {x,y) 
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X\ 


%2 


Sq(x), k = 1 


S (x), k = 2 


0.0 


0.0 


0.9170 


8950 


1.0 


1.0 


-0.0488 


-0 0630 


2 


9 


— Q 4^ x 1 0~ 3 


-0 0085 


3.0 


3.0 


—4.41 x 10~ 4 


3 30 x 10~ 3 


4.0 


4.0 


2 33 x 10~ Y 


1 56 x 1 0~ 3 


5.0 


5.0 


1 fi3 x 10~ 4 


5 96 x 1 0~ 4 


6 


fi 


1.12 x 10 -4 


— 2 52 x 1 0~ 4 


7.0 


7.0 


3.47 x 10~ 5 


— 1.72 x 10~ 4 


8.0 


8.0 


0.0000 


0.0000 


9.0 


9.0 


3.47 x lO" 5 


-1.72 x 10~ 4 


10.0 


10.0 


1.12 x 10" 4 


-2.52 x 10" 4 


11.0 


11.0 


1.63 x 10~ 4 


5.96 x 10~ 4 


12.0 


12.0 


2.33 x 10-'' 


1.56 x lO" 3 


13.0 


13.0 


-4.41 x 10~ 4 


3.30 x 10- 3 


14.0 


14.0 


-9.45 x 10~ 3 


-0.0085 


15.0 


15.0 


-0.0488 


-0.0630 



Table 2: The scalar part Sq(x) in the free fermion propagator of the higher- 
order Overlap Dirac operators, along the diagonal of a 16 x 16 lattice with 
ant iper iodic boundary conditions. 

Now we turn on a background U(l) gauge field with parameters ( hi = 0.1, 
h 2 = 0.2, 4 0) = 0.3, A { 2 0) = 0.4 and m = n 2 = 1, as defined in Eqs. (7) and (8) 
in Ref. ). Then we examine the behaviors of So(x,y) and S^(x, y) for the 
higher-order ( k = 1,2 ) Overlap Dirac operators. We find that both So(x,y) 
and 5*5 (x, y) are local in the higher-order ( k = 1,2 ) fermion propagators. In 
Table 3, we list the real parts and imaginary parts of S (x, 0) and S 5 (x, 0) for 
the second order ( k = 2 ) fermion propagator, along the diagonal ( x\ = x 2 ) 
of the 16 x 16 lattice. 

In general, the scalar part So(x,y) and the pseudoscalar part Ss(x,y) in 
the higher order ( k — 1,2 ) fermion propagators seem to be local, especially 
for near continuum gauge configurations. However, one cannot exclude the 
possibility that they may cause the perturbative instability of the pole of the 
fermion propagator |J. On the other hand, for the Neuberger-Dirac operator, 
we are sure that S (x,y) = a/25 X)y and S 5 (x,y) = for any gauge configura- 
tion, as well as the perturbative stability of the pole of the fermion propagator 
H]. So, from this viewpoint, the chiral properties of Neuberger-Dirac operator 
are better than those of higher-order Overlap Dirac operators. 
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Xl 


X2 


Re(So) 


Im(So) 


Re(S 5 ) 


Im(5 5 ) 


0.0 


0.0 


0.8950 


0.0 


2.94 x 10~ 3 


0.00 


1.0 


1.0 


-0.0609 


0.0156 


1.42 x 10~ 3 


-3.81 x 10~ 4 


2.0 


2.0 


-6.31 x 10~ 3 


5.97 x 10" 3 


-2.25 x 10~ 4 


1.30 x 10~ 4 


3.0 


3.0 


7.24 x 10~ 4 


-3.08 x 10~ 3 


-1.08 x 10~ 4 


3.48 x 10~ 4 


4.0 


4.0 


-8.25 x 10~ 4 


-1.43 x 10" 3 


6.01 x 10" 5 


9.28 x 10~ 5 


5.0 


5.0 


-6.54 x 10~ 4 


7.07 x 10~ b 


1.04 x 10" b 


-3.82 x 10~ b 


6.0 


6.0 


3.84 x 10~ 4 


-1.05 x 10~ 4 


-1.42 x 10~ b 


-1.30 x 10~ 5 


7.0 


7.0 


9.54 x 10~ 5 


-1.33 x 10~ 4 


8.16 x 10~ 6 


0.0 


8.0 


8.0 


2.14 x 10~ 5 


1.28 x 10~ 4 


0.0 


0.0 


9.0 


9.0 


2.27 x 10~ 4 


1.36 x 10~ 4 


-6.21 x 10~ b 


-2.31 x 10~ b 


10.0 


10.0 


1.49 x 10~ 4 


2.75 x 10~ 4 


-4.21 x 10~ b ' 


-1.34 x 10~ b 


11.0 


11.0 


-2.05 x 10~ 4 


-7.14 x 10~ 4 


2.17 x 10" b 


5.83 x 10~ b 


12.0 


12.0 


7.16 x 10~ 4 


-1.53 x 10~ 3 


-4.10 x 10~ b 


9.00 x 10~ b 


13.0 


13.0 


2.55 x 10~ 3 


-1.85 x 10~ 3 


-3.08 x 10~ 4 


1.76 x 10~ 4 


14.0 


14.0 


-8.30 x 10~ 3 


2.46 x 10" 3 


-2.70 x 10~ 4 


1.27 x 10~ 5 


15.0 


15.0 


-0.0629 


8.62 x 10~ 4 


1.46 x 10~ 3 


-2.99 x 10~ b 



Table 3: The scalar part So(x, 0) and pseudoscalar part S$(x, 0) in the fermion 
propagator of the second order ( k — 2 ) Overlap Dirac operator, in a back- 
ground gauge field ( see text for the parameters ) on a 16 x 16 lattice with 
antiperiodic boundary conditions. One of the end points of the propagator is 
fixed at the origin, while the other end point is at one of the sites along the 
diagonal ( x 1 — x 2 ). 
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4.2 Axial anomaly 

The axial anomaly of GW Dirac operator D satisfying ([1]) is [10, 1 1 

A L (x) = a tr[j 5 (RD)(x,x)} (70) 

where the trace runs over the Dirac and color space. Substituting R = 
(aj 5 D) 2k into ([FU|), we obtain 

A L (x)=tr[(a l5 D) 2k+1 (x,x)} . (71) 

The sum of the axial anomaly over all sites is equal to the index of D, 

J2A L (x)=n_-n + . (72) 

X 

If the index of D is equal to the topological charge Q of the gauge background, 
then the sum of the axial anomaly is equal to Q. However, it does not nec- 
essarily imply that Al(x) would agree with the topological charge density at 
each site. This happens only when D is local. 

Since the higher-order ( k > ) Overlap Dirac operator fl47|) is also topo- 
logically proper ( i.e., its index agrees with the background topological charge 
for any gauge background satisfying the topological bound ), then it follows 
that its axial anomaly would agree with the topological charge density at each 
site if D is local. ( i.e., the gauge configuration satisfies the locality bound 
which is more restrictive than the topological bound ). 

In the following, we compute the axial anomaly Al{x) in a two-dimensional 
background U(l) gauge field, for the Neuberger-Dirac operator and higher- 
order ( k — 1,2 ) Overlap Dirac operators respectively. We compare them to 
the topological charge density p(x) fl67|) of the gauge background on the torus. 

The deviation of the axial anomaly of a lattice Dirac operator in a gauge 
background can be measured in terms of 

1 ^ \A L (x) -a 2 p(x)\ 

where N s is the total number of sites of the lattice, and p(x) is the topological 
charge density inside the square of area a 2 centered at x. 

In a nontrivial U(l) gauge background with parameters Q = 1, hi = 0.1, 
h 2 = 0.2, 4 0) = 0.3, 4 0) = 0.4 and m = n 2 = 1 ( as defined in Eqs. (7) 
and (8) in Ref. [0 ) on the 12 x 12 lattice with a = 1, the axial anomaly 
Al(x) and its deviation 6 are computed for the Neuberger-Dirac operator and 
higher-order ( k = 1, 2 ) Overlap Dirac operators respectively. The results are 

0.110, Neuberger , 

0.193, k = l, (74) 
0.351, k = 2. 
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The relatively large deviations of axial anomaly in higher-order ( k = 
1,2 ) Overlap Dirac operators indicate that they are less localized than the 
Neuberger-Dirac operator. And the locality of D gets worse as the order k 
goes higher ( i.e., Dk =2 is less localized than D^ = i ). We have confirmed this 
by examining \D(x, y)\ versus \x — y\ explicitly. This implies that the e in the 
locality bound 

||1 — C/(p)|| < e, for all plaquettes (75) 

for a higher order ( k > ) Overlap Dirac operator is more restrictive ( smaller 
) than that of the Neuberger-Dirac operator, and it gets smaller as the order 
goes higher. 

For all background gauge configurations we have tested, the Neuberger- 
Dirac operator always gives anomaly deviation ( 5 ) smaller than those of the 
higher-order ( k = 1,2 ) Overlap Dirac operators. 



4.3 Fermion determinant 

The fermion determinant det(D) is proportional to the exponentiation of the 
one-loop effective action which is the summation of any number of external 
sources interacting with one internal fermion loop. It is one of the most crucial 
quantities to be examined in any lattice fermion formulations. The determi- 
nant of D is the product of all its eigenvalues 

det(D) = (1 + e k ) (n + +n - ) det'( J D) (76) 

where det'(D) is equal to the product of all non-zero eigenvalues. Since the 
eigenvalues of D are either real of come in complex conjugate pairs, det'(D) 
must be real and positive. For Q = 0, then n + + n_ = and det(D) = det'(D). 
For Q 7^ 0, then n + + n„ ^ and det(.D) = 0, but det'(-D) still provides 
important information about the spectrum. In continuum, exact solutions of 
fermion determinants in the general background U(l) gauge fields on a torus 
( L\ x L,2 ) was obtained in Ref. [12]. In the following we compute det (D) 



for the Neuberger-Dirac operator and higher-order ( k — 1,2 ) Overlap Dirac 
operators respectively, and then compare them with the exact solutions in 
continuum. For simplicity, we turn off the harmonic part ( h± — hi — ) 
and the local sinusoidal fluctuations ( Ai = A§ = ) and examine the 
change of det'(D) with respect to the topological charge Q. For such gauge 
configurations, the exact solution [|1^ 



is 



det'[D(Q)] = N ] j{^f l (77) 



where the normalization constant iV is fixed by 
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exact 


Neuberger 


k = 1 


k = 2 


Q 


det' [D(Q)]exact 


det'LD(Q)] 


det'fZXQ)U-i 

L \ / J fy — - 1 - 


det'[D(Q)k-2 

L \tS/J ft- — ^ 


i 


1 OOOOO 


1 00000 


1 00000 

J- . \J\J\J\J\J 


1 00000 


2 


2 82843 


2.77348 


3 01695 


3 05686 


3 


6.15840 


5.96891 


7.21380 


7.51909 


4 


11.3137 


10.7157 


14.6911 


15.7450 


5 


18.3179 


16.9340 


26.1688 


28.2029 


6 


26.8177 


24.2001 


41.9768 


44.4919 


7 


36.1083 


32.0006 


62.1673 


63.6961 


8 


45.2548 


40.2920 


86.8790 


85.2556 


9 


53.2732 


45.7353 


112.877 


106.173 


10 


59.3164 


50.2816 


141.546 


128.122 



Table 4: The fermion determinant versus the topological charge Q. The nor- 
malization constant is chosen such that det'[D(l)] = 1. The results on the 8x8 
lattice are listed in the last three columns for the Neuberger-Dirac operator 
and higher-order ( k — 1, 2 ) Overlap Dirac operators respectively. The exact 
solutions on the 8x8 torus are computed according to Eq. (ff7|). 



such that det'[D(l)] = 1. 

In Table 4 and Table 5, the fermion determinants det'(D) are listed for 8 x 8 
and 16 x 16 lattice respectively. It is clear that the Neuberge-Dirac operator 
always produces results better than those of the higher-order ( k = 1,2 ) 
Overlap Dirac operators. The first order ( k = 1 ) Overlap Dirac operator 
performs better than the second order ( k — 2 ) one. This is essentially due to 
the fact that D becomes less localized as the order ( k ) goes higher. 



5 Discussions and Conclusions 

We can understand the emergence of Fujikawa's proposal by the following 
considerations. 

If one requires that D is 75-hermitian, 

D ] = 75^75 , 

and normal, 

D ] D = DD ] , 

( Note that these two conditions are sufficient to guarantee that the real eigen- 
modes of D have definite chirality (|28|) and satisfy the chirality sum rule ([31]), 
and each complex eigenmode has zero chirality (27). ), then one immediately 
obtains 

75.075.0 = .D75-D75 . 
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exact 


Neuberger 


k = 1 


k = 2 


Q 


det' [D(Q)]exact 


det'LD(Q)] 


det'fZXQ)U-i 


det'[D(<9)lfe- 2 

L \tS/J ft- — « 


i 


1 OOOOO 


1 00000 

J- . \J\J\J\J\J 


1 00000 

J- . \J\J\J\J\J 


1 OOOOO 


2 


5 65685 


5 66186 


5.76472 


5 78270 


3 


24.6336 


24.0615 


25.8620 


26.1218 


4 


90.5097 


90.4894 


98.9598 


100.627 


5 


293.086 


286.003 


337.077 


346.382 


6 


858.166 


822.664 


1046.70 


1092.03 


7 


2310.93 


2170.94 


3018.17 


3192.52 


8 


5792.62 


5354.28 


8198.43 


8786.52 


9 


13637.9 


12309.5 


20911.5 


23111.2 


10 


30370.0 


27336.6 


51335.5 


57563.2 



Table 5: The fermion determinant versus the topological charge Q on the 
16 x 16 lattice. The normalization constant is chosen such that det'[Z)(l)] = 1. 
The results on the lattice are listed in the last three columns for the Neuberger- 
Dirac operator and higher-order (k = 1, 2) Overlap Dirac operators respec- 
tively. The exact solutions on the 16 x 16 torus are computed according to 
Eq. (0). 

Multiplying above equation by 75, we obtain 

= 75^75^75 , 

which can be rewritten as 

7 5 (a7 5J D) 2 = (a7 5J D) 2 7 5 . 

Since (a'y^D) 2 is Hermitian and commutes with 75, one finds an example of 
R = (a'y^D) 2 which depends on D. Then it is straightforward to generalize 
this R to any powers, 

R = (a l5 D) 2k , fc = 0,l,2,... 

Substituting this R into the GW relation ([!]), we obtain (|3]) which is equivalent 
to Fujikawa's proposal (Q). 

It seems to us that Fujikawa's higher-order realization of the Overlap Dirac 
operator may not be feasible for practical computations in lattice QCD, in view 
of its locality, chiral properties and computational accessibility in comparison 
with those of the Neuberger-Dirac operator. Nevertheless, from a theoretical 
viewpoint, it has widened our scope and deepened our understanding of the 
Overlap which does capture one of the fundamental aspects of the nature. 
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